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ABSTRACT 



We derive analytical solutions and dispersion relations of global magnetic 
Poincare (magneto-gravity) and magnetic Rossby waves in the approximation of 
shallow water magnetohydrodynamics. The solutions are obtained in a rotating 
spherical coordinate system for strongly and weakly stable stratification sepa- 
rately in the presence of toroidal magnetic field. In both cases magnetic Rossby 
waves split into fast and slow magnetic Rossby modes. In the case of strongly 
stable stratification (valid in the radiative part of the tachocline) all waves are 
slightly affected by the layer thickness and the toroidal magnetic field, while 
in the case of weakly stable stratification (valid in the upper overshoot layer 
of the tachocline) magnetic Poincare and magnetic Rossby waves are found to 
be concentrated near the solar equator, leading to equatorially trapped waves. 
The frequencies of all waves are smaller in the upper weakly stable stratification 
region than in the lower strongly stable stratification one. 

Subject headings: Sun: interior - Sun: magnetic fields - Sun: oscillations 
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Introduction 



The tachocline, w hich is believed to exist near the base of the solar convection zone 
1^ 



( ISpiegel &: Zah 



last years (iGough 



9921). has become a subject of rapidly growing scientific interest in the 



20071 ). The existence of the tachocline is important for two main reasons. 



First, it may preven t the spread of the so 



zone to the interior (ISpiegel fc Zahn 



1992 



ar angular momentum from the convec tion 



Gough fc Mclntvr 



1998 



Gar and 



20071 ). thus 



retaining the inferred solid body rotation in the radiative zone; second, the processes of solar 
magnetic field generation and global magnetohy drodynamic (MHO) instabilities, which ar e 



cruci al for solar activity, may take place there (jPikpati fc Gilman 



2005 



Gilman fc Gaily 



20071 ). MHD waves and oscillations in the tachocline may play a significant role in both 
cases: they may redistribute the angular momentum in the horizontal direction and some 
of them may become unstable due to the differential rotation, leading to magnetic flux 
emergence at the solar surface. Therefore, to study the complete spectrum of possible wave 
modes in this system is of vital importance. The tachocline is very thin compared to the 
solar radius, therefore the ordinary shallow water approximati on modified by the presence 



of a horizontal large-scale magnetic field can be easily applied (j Gilman 



2000). 



The spectrum of various shorte r scale shallow water MHD waves has been recently 



studied in Cartesian coordinates by 



Schecter et al. 



(1200 ll ). However, global wave modes. 



those with a waveleng th comparable to t h e sola r radius, must be considered in spherical 



coordinates. Recently, 



Zaqarashvili et al. 



(120071 ) have studied the spherical shallow water 



MHD waves in the simplest case, in which the high frequency branch i.e. magnetic Poincare 
waves (or magneto-gravity waves) and the influence of the tachocline thickness on the wave 
dynamics have been ignored (that is, when the surface gravity speed is much higher than 
the surface rotation speed). This approximation has enabled us to study the dynamics of 
global magnetic Rossby waves in the lower, strongly stable part of the tachocline, but it 
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fails in the upper, weakly stable overshoot region, where a ne gative buoyan cy due to the 



subadiabatic stratification strongly reduces the gravity speed (iGilman 



2OOOI) 



In this letter, we present the analytical spectrum of global linear shallow water MHD 
waves for both parts of the tachocline in the absence of differential rotation. 



2. Global shallow water MHD waves 



We use the linearized shallow water MHD equations in the rotating spherical coordinate 
system (r,9,(t)), where r is the radial coordinate, 9 is the co-latitude, and is the longitude 



(IZaqarashvili et al. 



20071 ). We also assume the unperturbed toroidal magnetic field to be 



= B(9) sin^, where B{9) may have an arbitrary profile with latitude, which can be 
chosen later. Then a nor mal mode analysis of the form exp(— zcjt + iscj)) gives the equation 



(IZaqarashvili et al 



20071 ) 



{\D + s^) 



s2-eA2(l-/i2) 



a 



a 



+s a ue + 2a fiDiiQ + fis—{D + 2jj,)u0 + j^D{a )ue = 0, 



A 



(1) 



with 



A 



UJ 



2VLn 



va 



/i = cos 6', D 



d_ 



d 



sin 6'—, Ue = sinOug, 
o9 



where ug is the poloidal velocity component, g is the reduced gravity in the tachocline, p 
is the medium density, Qq is the angular velocity, Rq is the distance from the solar center 
to the tachocline, and Hq is the tachocline thickness. This is the general equation for 
the linear dynamics of the shallow water MHD system, which contains various kinds of 
waves (magnetic Poincare and magnetic Rossby waves). When e = and B = const, this 
equation transforms into the associated Legendre equation and governs the dynamics of fast 
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and slow magnetic Rossby waves (IZaqarashvili et al.l 120071 ). This approximation neglects 
magnetic Poincare waves and the influence of the tachocline thickness on the dynamics of 
magnetic Rossby waves. Therefore, equation (1) must be solved with e 7^ in order to give 
the complete spectrum of waves. The subadiabatic stratification in the tachocline provides 
negative buoyancy, which leads to a reduced gravity and consequently an increased value 
of e. In the lower strongly stable part of the tachocline e still is much smaller than unity, 
but in upper overshoot layer e becomes much greater than unity. Therefore, the solution of 
equation (1) for the two extreme cases e ^ 1 and e ^ 1 covers both parts of the tachocline. 
Here we use a magnetic field profile typical of the Sun, namely 



(2) 



19971 ). 



which has opposite signs in the northern and southern hemispheres ( iGilman fc Foxl 
The magnetic field profile (2) leads to a = aofi, where ao = Bq/ {2VIqRq^^ At: p) . In the next 
subsections we will give analytical solutions and dispersion relations of global shallow water 
MHD waves for each case separately. 



2.1. e ^ 1 (strongly stable stratification; valid for the radiative pcirt of 

tachocline) 



Reduced gravity in t 
as 500-1.5 ■ 10"^ cm-s~^ ( ISchecter et al 



le strongly stable ra diative part of the tachocline can be estimated 
2OOII ). Then, for the layer thickness HQ=l{f cm 



we get e = 4.5 -10 - 0.13 (where ^0=2. 6 -10 ® s ^, Ro=5- 10^° cm have been used). 
Therefore, e <^ 1 is a good approximation in the radiative part of the tachocline. 

In this approximation and using the weak magnetic field limit, i.e. <^ 1, equation 



(1) leads to the spheroidal wave equation (second order terms with e and an are neg. 



whose typical solutions are the spheroidal wave functions Ssn (lAbramowitz fc Stegun 



ected). 



mi), 
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where n plays the role of the poloidal wavenumber. The approximate solution can be 
written as 



/ .^-eA^(l-/i^) 

Ug OC Uobsn{(^i,IJ)\l — ^ ^ ^ , 



(3) 



where Uq is the amplitude and 



Note that n is integer for the solutions, which vanish at the poles (/i = ±1)- The dispersion 
relation is found to be 



A 



n[n + 1) 



A 



+ = 0. 



(4) 



The dispersion relation can be split into high (magnetic Poincare) and low (magnetic 
Rossby) frequency branches. For the magnetic Poincare waves we have 

AmP 



n(n + 1) 



11 



(5) 



This is the dispersion relation of ordinary su rface gravity waves and we can recover it for 



the high s harmonics (jLonguet-Higgins 



For the magnetic Rossby waves we have approximately 

2 



n{n + 1) 



A 



A 

s 



a 







0. 



(6) 



It is clear from equation (6) that 



;he magnetic field cause s the splitting of ordinary Rossby 



waves into fast and slow modes (jZaqarashvili et al. 



20071 ). For ao -C 1 the dispersion 



relations for the fast and slow magnetic Rossby modes can be approximated as 



A 



/ 



A, 



a, 



0' 



(7) 



s n{n + 1) " s 
We see that the fast magnetic Rossby mode has a dispersion relation similar to that of the 
ordinary Rossby waves but slightly modified by the magnetic field. Slow magnetic Rossby 
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waves, however, are new wave modes as their dispersion relation is different from Rossby 
and Alfven wave dispersion relations (the dispersion relation of Alfven waves should be 
= ±sao). There is a significant difference between the frequencies of fast and slow 
magnetic Rossby modes for sufficiently small ao- 

The dispersion diagram for the s = 1 harmonics of shallow water MHD waves 
according to Eq. (jl]) is shown in Fig 1. Here e = 0.01 and ao = 0.05. The frequencies of 
magnetic Poincare waves (solid lines) are much higher than the rotational frequency, Qq, 
and they increase with increasing n. For example, the frequency of the n = 3 harmonics 
is ~ 50^0) which corresponds to oscillations with period ~ 13 hours. On the other hand, 
the frequencies of magnetic Rossby waves are in general much lower than the rotational 
frequency. The absolute value of the frequency of fast magnetic Rossby waves (dotted line 
in the zoom) decreases with increasing n, such as happens with the HD Rossby wave. 

Fig. 2 shows the dependence of the s = l,n = 2 and s = l,n = 3 harmonics on e and 
ao (i.e. the magnetic field strength). The frequency (in absolute value) of fast (dotted) 
and slow (dashed) magnetic Rossby modes significantly increases when the magnetic field 
is increased. We also see that the frequency difference between the n = 2, 3 harmonics of 
fast magnetic Rossby waves is almost independent of ao. On the contrary, the n = 2, 3 
harmonics of slow magnetic Rossby waves diverge when the magnetic field is increased. It 
turns out that the frequencies of magnetic Rossby waves have almost no dependence on e, 
therefore they are not shown in the top panel. On the other hand, the frequency of magnetic 
Poincare waves depends on e (lower panel) and significantly decreases with increasing e. 
However, these frequencies have almost no dependence on oq, therefore magnetic Poincare 
waves are only slightly affected by the magnetic field. 
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2.2. e ^ 1 (weakly stable stratification; valid in the upper overshoot part of 

tachocline) 

In the weakly stable over shoot region (or say u pper tachocline) reduced gravity is much 



200l|). Then, for the thickness Ho=5 ■ 10^ cm 



smaller, being 0.05-5 cm-s ^ (iSchecter et al. 
we get e = 27 - 2.7 • 10^. Therefore, e ^ 1 is a good approximation in the overshoot region. 



In this case, we follow the calculation of iLonguet-Higgind (119681 ) (this paper considers 
spherical hydrodynamic shallow water waves in the Earth context), and introduce a new 
function i] = e^/'^jj, and also consider that A = e^^/^L and ao = e^^/^L, where L is of order 
unity. Then after some algebra and keeping only large terms equation (1) leads to the 
Weber (parabolic cylinder) equation, whose solutions finite a,t t] ±cx3 can be expressed in 
terms of Hermite polynomials Hy{ri) of order z/, where v satisfies 



2„2 



eA^ 



A A2 

and the solution is proportional to 



(2z/ + l)v^, 



(z/ = 0,l,2,...) 



1 2 

U0 cx e~2'' Hy{ri). 



(9) 



Therefore the function ug is exponentially small beyond the turning points cos 6* = 
g-i/4^2zy + 1. Due to the large e, the co-latitude 6 satisfying this condition is close to 90°, 
therefore the solution is confined to the neighborhood of the equator. 



Equation (8) yields the wave dispersion relation as 



A^-^A^-^A + 



2^2 



0. 



v^e e e 

n the non-magnetic case this equation transforms into equation (8.11) of 



(10) 



Longuet-Higgins 



f 19681 ). This is a fourth order equation and it describes magnetic Poincare and magnetic 
Rossby waves. For large values of e, the solutions can be written separately for both types 
of waves. 
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The solutions for Poincare waves can be written approximately as 
It turns out that the magnetic field has almost no infiuence on the dynar nics of Poincare 



wave s and the dispersion relation is the same as for the nonmagnetic case (jLonguet-Higgins 



19681). 

The dispersion relation for the magnetic Rossby waves can be written as 

{2u+l)^/^\^ + s\- als'^ = 0. (12) 

The magnetic field causes the splitting of ordinary Rossby waves into fast and slow magnetic 
Rossby waves as in the e <^ 1 case. If y/ea^ <^ 1, then we can write 

Xf 1 9 Ao 9 

Fast magnetic Rossby waves are similar to hydrodynamic Rossby waves modified by the 
magnetic field. Slow magnetic Rossby waves have a similar dispersion relation as in the 
e ^ 1 case. 

The dispersion diagram for the s = 1 harmonics of magnetic Poincare and magnetic 
Rossby waves according to Eq. (10) is shown in Fig 3. Here e = 2700 and ao = 0.05. The 
behaviour of fast magnetic Rossby waves (dotted line) is similar to the small e case: the 
absolute value of their frequency significantly decreases with increasing u (which now plays 
the role of the poloidal wave number). However, we immediately note that the frequency 
now is much smaller than that of the e -C 1 case. This is also true for magnetic Poincare 
waves (solid lines): their frequency is significantly small compared to the case of small 
e, now being even smaller than the rotational frequency. The frequency of slow magnetic 
Rossby waves only slightly depends on z/, but now it decreases with increasing i>, contrary 
to e ^ 1 case. 
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Fig. 4 shows the dependence of the s = l^u = 2 and s = l,u = 3 harmonics on e 
and ao- The frequency of magnetic Rossby mode harmonics significantly depends on ao; 
the frequencies of both fast and slow magnetic Rossby modes increase with increasing ag. 
On the other hand, it turns out that the frequencies of magnetic Poincare waves have 
almost no dependence on ao- However, they are significantly reduced with increasing e, 
becoming smaller than the rotational frequency for e > 1000. Fast magnetic Rossby waves 
also depend on e, their frequency significantly decreasing when this parameter is increased. 
Slow magnetic Rossby waves have almost no dependence on e as it is suggested from Eq. 
( |T3l) . Thus increasing e, which is equivalent to reducing g, generally leads to a decrease of 
the wave frequencies. 

The numerical solution of the dispersion relation (10) shows that the frequency of the 
magnetic Rossby wave harmonic with u = becomes complex for sufficiently large i-e- 
large magnetic field. This may point to some kind of instabilities similar to the polar kink 
instability found by Gaily (2003). However, the value of ao for which the frequency becomes 
complex is significantly outside the range for which Eq. (10) is valid (note that ao should 
be proportional to e~^/^). Therefore, the complex solution can be spurious and may not 
refiect a real instability. 

3. Conclusions 

We have derived analytical solutions and dispersion relations of global shallow water 
MHD waves for the solar tachocline in a rotating spherical coordinate system. The 
solutions and dispersion relations are obtained for weakly and strongly stable stratifications 
separately. The weakly stable stratification is valid in the upper overshoot part of tachocline, 
while the strongly stable stratification is valid for the radiative part of tachocline. The 
solutions include the Poincare and Rossby waves modified by the presence of the magnetic 
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field. 



Here we use a realistic latitudinal profile of the toroidal magnetic field ~ cos 6 sin 9 



that changes sign at the equator. Similar to the e = case (jZaqarashvili et al. 



20071 ). the 



magnetic field leads to the splitting of hydrodynamic Rossby waves into fast and slow 
magnetic Rossby modes. 



In the case of strongly stable stratification, the 



waves is similar to that of the hydr odynamic case (ILonguet-Higgins 



relation of magnetic Rossby waves (jZaqarashvili et al 



dispersion relation of m agnetic Poincare 



19651 ). The dispersion 



20071 ) is slightly modified by the 



thickness of the layer. In the case of weakly stable stratification, both magnetic Poincare 
and magnetic Rossby waves are concentrated near the equator, thus leading to equatorially 
trapped waves. 

Another remarkable feature is that the frequencies of all waves are significantly reduced 
in the large e case compared to the e <^ 1 one. This means that the upper overshoot part 
of the tachocline supports oscillations with smaller frequency than the lower radiative part 
of the tachocline. 

The solutions have been obtained without taking into account the differential rotation 
which is present in the tachocline. The future inclusion of differential rotation probably will 
l ead to the instability of some modes, which was found by simulations in the inertial frame 



flGilman fc Callv 



20071 . and references therein). 
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Fig. 1. — Dispersion diagram of spherical shallow water waves in the presence of a toroidal 
magnetic field for small e. The two extreme upper and lower solutions (solid lines) correspond 
to magnetic Poincare waves, whereas the two low frequency modes shown also in the zoom are 
fast (dotted) and slow (dashed) magnetic Rossby waves. Symbols denote solutions obtained 
with the approximate formulas ([5]) and ([7]). The parameters used to obtain the dispersion 
diagram are e = 0.01, s = 1, and ao = 0.05. 
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Fig. 2. — Upper panel: frequencies of the s — l,n — 2 and s — l,n — 3 harmonics of 
fast (dotted) and slow (dashed) magnetic Rossby modes vs. ao for e = 0.0045. Magnetic 
Poincare waves do not depend significantly on the magnetic field, therefore they are not 
shown here. Lower panel: frequencies of the s — l,n — 2 and s — l,n — 3 harmonics of 
magnetic Poincare waves vs. e for ao — 0.05. 
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Fig. 3. — Dispersion diagram of spherical shallow water waves in the presence of a toroidal 
magnetic field for large e. The two extreme upper and lower solutions (solid lines) correspond 
to magnetic Poincare waves, whereas the low frequency modes (also shown with zoom) are the 
fast (dotted) and slow (dashed) magnetic Rossby waves. Symbols denote solutions obtained 
with the approximate formulas ffTTl) and f|T3l) . The parameters used to obtain the dispersion 
diagram are e = 2700, s = 1, and ao = 0.05. 
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Fig. 4. — Upper panel: frequencies of the s — l.,v — 2 and s — — ?> harmonics of 
fast (dotted) and slow (dashed) magnetic Rossby modes vs. ckq for e = 400. Middle panel: 
frequencies of the s = 1, z/ = 2 and s = l,z/ = 3 harmonics of magnetic Poincare and 
magnetic Rossby modes vs. e for ao — 0.05. Lower panel: same as in middle panel, but only 
for magnetic Rossby waves. 



